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Abstract 

It is a common belief among field theorists that path integrals can be com- 
puted exactly only in a limited number of special cases, and that most of 
these cases are already known. However recent developments, which gener- 
alize the WKBJ method using equivariant cohomology, appear to contradict 
this folk wisdom. At the formal level, equivariant localization would seem 
to allow exact computation of phase space path integrals for an arbitrary 
partition function! To see how, and if, these methods really work in prac- 
tice, we have applied them in explicit quantum mechanics examples. We 
show that the path integral for the 1-d hydrogen atom, which is not WKBJ 
exact, is localizable and computable using the more general formalism. We 
hnd however considerable ambiguities in this approach, which we can only 
partially resolve. In addition, we find a large class of quantum mechanics 
examples where the localization procedure breaks down completely. 



1. Introduction 

The major failing of the path integral method applied to quantum field theory or 
quantum mechanics is that most physically relevant path integrals cannot be com- 
puted. The known exceptions fall roughly into three categories. There are the Gaus- 
sian forms, which are the basis of perturbation theory in quantum field theory. There 
are also certain exactly solvable models, usually quantum mechanical integrable sys- 
tems. In such cases the path integral is often computable due to the fact that the 
WKBJ approximation is exact. Thus the path integral localizes to a sum over clas- 
sical trajectories. The third category is path integral representations of topological 
quantities such as the Euler character or the Witten index. Here one often finds that 
supersymmetry in the path integral implies an equivariant cohomology structure, 
which localizes the path integral onto zero modes or characteristic classes. 

As emphasized in [jl], || |3], the common denominator of most exactly computable 
path integrals is equivariant cohomology and localization. This raises the possibil- 
ity that by studying these methods directly we may greatly expand our inventory of 
doable path integrals, as well as WKBJ-like approximation schemes. Though equiv- 
ariant cohomology is usually associated with supersymmetric theories, it actually has 
very general application. The key ingredient required is some sort of differential form 
structure among the physical, auxiliary, or ghost variables. This is a rather weak 
requirement. In fact, as recently shown by Blau, Keski-Vakkuri, and Niemi]I|], and 
by Niemi and TirkkonenQ], equivariant cohomology can be exhibited in an arbitrary 
phase space path integral. Further, they showed formally that under seemingly weak 
conditions, this results in equivariant localization of these path integrals. They would 
thus be exactly computable! 

Of course, we do not really expect arbitrary phase space path integrals to be 
exactly computable. Since equivariant localization proofs are very formal, there are 
ample opportunities for subtle difficulties to creep in. Also phase space path integrals 
are themselves rather disreputable. As pointed out in [||], there is a well known 
false method for computing arbitrary phase space path integrals via Hamilton- Jacobi 
theory and (invalid) canonical transformations of the path integral variables |5|. |(|. 

In short, there is no substitute for simple, explicit examples, and that is the 
subject of this paper. In section 2. we briefly review equivariant cohomology and 
localization formulae for path integrals. In section 3. we make the formalism explicit 
by considering the phase space path integrals for the free particle and the harmonic 
oscillator. We use equivariant localization to write the partition function of a one- 
dimensional quantum mechanical system with an arbitrary static potential as an 
elementary contour integral. Along the way we see that this derivation breaks down 
for essentially all potentials which are bounded below. 

In section 4. we consider the 1-D hydrogen atom, whose potential is not bounded 
below. Its path integral is not WKBJ exact, and is not amenable to Morse theory 
analysis since its classical trajectories all coalesce at two points in phase space. Nev- 
ertheless we show that it is exactly computable via equivariant localization, modulo 
an important ambiguity which we discuss. 
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2. Equivariant Cohomology and Localization 



Equivariant cohomology is a simple and powerful generalization of ordinary de 
Rham cohomology and the calculus of differential forms (for details, see 0, ||). To 
describe differential forms with respect to some rf-dimensional manifold M, it is con- 
venient to introduce a contravariant vector of anticommuting Grassmann variables 
c a . Then differential forms can be represented as covariant tensor functions on M 
contracted with the c a 's: 



— forms 

1 — forms 

2 — forms 



<Pa{x)c a 

<p ab (x)c a c'\ 

etc.. 



The exterior derivative operator d can be written 

d ="h (2 - 2) 

It takes Ti-forms to n+l-forms and is nilpotent: d 2 = 0. De Rahm cohomology is the 
classification of closed forms (forms which are annihilated by d) modulo exact forms 
(which are d of something). 

To define equivariant cohomology, we first introduce a new operation, i x , which is 
interior multiplication with respect to an (arbitrary) vector x a { x )- Explicitly: 



i x = J dc aX a (2.3) 
Clearly % x takes an n-form to an n— 1 form: 

i x (f>(x) = 0, 

i x $ a {x)c a = <p aX a , (2.4) 
t x cj )ab (x)c a c b = ( p ab (x a c b -x b c a )- 

Furthermore % x is nilpotent. 

We can now define an equivariant exterior derivative d x by 



d x = d + i x . (2.5) 
m n-form to both 

generally nilpotent: 



Note that d x takes an n-form to both an n+l-form and an n—l form. Now, d x is not 



d x = (di x + i x d) = £ x , (2.6) 



where C x is the Lie derivative with respect to the vector x a { x )- Thus d x is nilpotent 
only on the subspace of differential forms which are annihilated by the Lie derivative 



C x . 
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Without going into more rigorous mathematical detail, we are already in a position 
to sketch the general derivation of equivariant localization formulae. Suppose we want 
to compute an integral over the manifold M of some general form a, where by a general 
form we mean a linear combination of different n- forms (including a top-form): 



M 



a (2.7) 



In order to apply equivariant localization, we must be able to find a vector x a such 
that, with respect to x a , ol is equivariantly closed: 

d x a = . (2.8) 

In addition, we must be able to find some other differential form (3 which is in the 
subspace annihilated by C x . C x (3 = 0, and is not equivariantly exact. When this can 
be done, one constructs the following modified integral: 



a e~ xd xP , (2.9) 



M 

where A is a numerical parameter. It is then straightforward to show that, formally, 
this modified integral is independent of the value of A: 

± / a e~ M ^ = -A / (dJ3)a e~ Xd ^ 
d\JM Jm x 

= A 



M 



. 



In the above we used d x a=0, d x d x /3=0, and the fact that the integral of an equivari- 
antly exact form vanishes. 

Thus, provided the limits A^O and A^oo exist, we obtain the following localiza- 
tion formula: 

-\d x f3 



a= lim / ae- Xdxli . (2.11) 

M A^oo Jm 

The original integral over M has been localized to a subspace M x , which is the support 
for the nontrivial equivariant cohomology {d x (3 = 0). 

When applied to ordinary (finite dimensional) integrals over compact manifolds, 
the above arguments can be made rigorous, and the right hand side of (|2.11|) (for 
suitable choice of (3) reduces to the WKBJ formula. This is the Duistermaat-Heckman 
theorem. 

At the formal level, this result was extended in |3J to apply to general phase space 
path integrals. Consider a general bosonic quantum mechanical phase space path 
integral: 

Z{T) = J [d(Liouville)] e iS = J [d<p a ] y / detjjo^jexp % £ dt (ej a - H (</>)) , (2.12) 



4 



where <fi a (t) are the coordinates of phase space (satisfying periodic boundary condi- 
tions), 9 a {4>) are their conjugates, H is the hamiltonian, and u ab is an antisymmetric 
tensor giving the symplectic 2-form whose matrix inverse defines Poisson brackets: 

Uab = d a 6 b - d b 9 a . (2.13) 

Introducing real Grassmann auxiliary coordinates c a (t), which also satisfy periodic 
boundary conditions, we can write: 

Z(T) = J [d<f) a ] [dc a ] exp i J* dt (6 a <p a - H(</>) + \c a u ab c b ) . (2.14) 

We can then identify loop space differential forms as tensor functionals of the a (£)'s 
contracted with the c a (£)'s. There are two vectors of obvious significance for defining 
equivariant operations. One is the hamiltonian vector field, which is defined by the 
relation: 

d a H{4>) = u abX b H . (2.15) 

The other is <p a , since <p a d a generates time translations on our differential forms. 
In M, equivariant cohomology is defined with respect to the vector 

x a s (t) = rit) - xum , (2.i6) 

since, as one can easily verify, dsS = 0. Furthermore, the zeroes of correspond 
to the classical trajectories. Thus, if we shift the action by d$ of a 1-form f3 

S^S + XdsP, (2.17) 

and choose /3 proportional to x a si we can expect to obtain a WKBJ localization of 
Z(T). Since both Xs an< i ° a are contravariant, we cannot construct such a (3 without 
also introducing a metric on loop space, G ab (ti, £ 2 ). Then we may define (3 as 

(3 = J dhdt 2 G ab (t u t 2 )x a s {tx)c b {t 2 ) . (2.18) 

Recalling that the localization argument requires £s/3=0, one finds that this reduces 
to the condition 

C s G ab (t ll t 2 ) = 0. (2.19) 

In other words, Xs m ust be a Killing vector of the loop space metric. It is easy to see 
that automatically annihilates G ab (ti,t 2 ) provided it takes the form 

G ab (t u t 2 ) = g ab (0(t)) 8{t x - t 2 ) . (2.20) 

Then ( |2.19| ) reduces to the constraint 

Ch9M) = • (2.21) 
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In other words, Xh m ust be a Killing vector of the finite dimensional metric g a b- 

As shown in the formal procedure just outlined results in a WKBJ localization 
of the path integral ( 2.12 ). Thus in the case where these formal steps actually carry 
through, the path integral is WKBJ exact. As already mentioned, there are a number 
of known examples where this occurs. This is well-trodden territory, and we have not 
attempted to explore it further. 

Instead, we turn to the recent paper by Niemi and Tirkkonen[|j], in which they 
apply a different equivariant localization scheme to the general phase space integral 
Q2.12| ). The difference lies in the choice of the 1-form (3: instead of ( |2.18| ) they choose 



(2.22) 



where we suppressed the t integrations which accompany index contractions. The 
equivariant localization procedure is otherwise the same as described above. One 
then obtains a localization of Z(T) not onto classical trajectories, but rather onto the 
time independent modes of <p a (t) and c a (t). Their result (correcting a minor typo) is 



Z(T) 



%dCQ Q~ iTH ^ +i 2 Tc » UJabC 



\ 



Det 



ab) 



sinh 



{Qab + Rab) 



(2.23) 



In this expression, the determinant is an ordinary matrix determinant, and the inte- 
grals are finite dimensional integrals over the time-independent real variables 0g and 
real Grassmann variables eg. The tensors Q a b and R a b can be expressed as geometric 
quantities in terms of the metric g a b and the hamiltonian vector x°k which is a Killing 
vector of g ab : 



Rab — Rabcd,C()CQ 



Now ( |2.23j ) is a remarkably simple result. If this localization proved to be valid 
for some large class of physically relevant partition functions, it would be an excit- 
ing and important achievement. Less optimistically, we have here at the very least 
an alternative to WKBJ methods which is equally nonperturbative. Further, both 
WKBJ and the Niemi-Tirkkonen formula are merely examples of much more general 
equivariant localization techniques. 



3. Application to 1-D Quantum Mechanics with Static Potentials 

The simplest possible application of ( |2.23| ) is to one dimensional quantum mechan- 
ics with static potentials. In this case there are only two phase space coordinates, 
x(t) and p(t), and the Hamiltonian is of the form 



H 



1 - 2 + V(x). 



(3.1) 
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For the moment, let us confine ourselves to potentials which are bounded below. Then 
it is convenient to define "harmonic" and polar coordinates as follows: 



p = rsm9, 
H=\p 2 + \u?y 2 + V min = \r 2 + V m 



V(x) - V min = \u 2 y 2 = rcostf, (3.2) 



Note that we only change variables in the time independent coordinates of the local- 
ized integral, never in the original phase space path integral. 

The advantage of working in polar coordinates is that the hamiltonian vector field 
has only one nonvanishing component: 

x = x = - u Tx ' (3 ' 3) 

and it is straightforward to solve the constraint C,H9ab = 0, which is equivalent to: 

xdeQee + ^QeedeX = , 

de(x9re) + 9eed r X = , (3.4) 
Xdegrr + 2g r edrX = ■ 

The general solution is: 

_ fir) 
9ee — — — , 

f( r ) 

Note that the metric is not completely determined by the constraint; the general 
solution (3.5) involves three arbitrary functions of r: /, f 2 , and fa. This is not 
surprising, since we expect that the localized integral (|2.23|) is at least partially metric 
independent. 

A much worse feature of the general solution (3.5) is that g^ir, 0) is not, in general, 
a single- valued function. To see this, observe that g r g(9=0) = g r o(6=2'n) only obtains 
if 

d r - = , (3.6) 
Jo X 

which is equivalent to 

C 2-7T £jf<£ 

d6 — = constant . (3.7) 
o ay 



The harmonic oscillator gives the only solution of (|3.7|) . Thus we conclude that, 
for essentially all static potentials which are bounded below, the Niemi-Tirkkonen 
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equivariant localization procedure fails, due to the nonexistence of a single-valued 
metric satisfying the Lie derivative constraint. 

What of the harmonic oscillator, the sole survivor of this calamity? In this case 
we have 

X= -u > 

9ee = — t- , 9re = , fiW = -r + fa{r) . 

The antisymmetric tensor VL^ and the curvature tensor R a b c d each have only one 
independent component: 

O - //(r) R ~ c/ ( r ) n rqn\ 

cj 2wc(r) 

where we have introduced the (as yet) arbitrary function c(r): 

It is now straightforward to plug into the localization formula ( 12.231) . Performing the 
Grassmann integrals and the 9 integral leads to the following simple result: 

Z ( T) = ir4(^). (3.11) 
u Jo ar Vsm cl J 

So, in the special case of the harmonic oscillator, the localization formula reduces to 
a total derivative. The final result is 

Z(T) = - r ( °Li • (3-12) 
v ; wsin[c(0)T] v ; 

It is surprising that, even in such a simple the harmonic oscillator, the 

final result of the localization procedure is ambiguous, i.e. it depends on the value 
at r=0 of the arbitrary function c(r). The correct result for the path integral is only 
obtained if we impose the additional ad hoc boundary condition 

lim c(r) = \uj . (3.13) 

As we will see in the next section, such ambiguities appear to be a general feature of 
equivariant localization of phase space path integrals. Thus the localization formulae, 
even when they are otherwise valid, are incomplete. This could have been anticipated 
even in the formal derivation. An obvious requirement which must be incorporated 
into the localization procedure is that the 1-form (3 should be homotopically equivalent 
to under the "supersymmetry" transformation generated by ds- Thus one expects 
generally that additional inputs are required when choosing g a b, in order to ensure 
that we are in the trivial homotopy class. 
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Before moving on to more interesting problems, let us consider the phase space 
path integral for the free particle. In this case H = \p 2 ■, and x x= X = P- The localiza- 
tion procedure again has ambiguities, which can be eliminated by fixing the metric 
for the free particle to be: 

«*=(£ t:Mo ?)■ < 3 - i4 » 

The localization formula can be reduced to 

(3.i5) 



With the metric (|3.14j) we have c(p) vanishing identically, and the above produces 
the correct result. 

4. The 1-D Hydrogen Atom 



The 1-d hydrogen atom is the static potential problem with hamiltonian|TO 

c 2 

■n 2 — - 

2 



H = \P 2 - h ■ (4-1) 



x 



While the partition function can be computed exactly by solving the Schrodinger 
equation, the path integral for the partition function cannot be directly performed 
by standard techniques. In particular it is not WKBJ exact. The classical bound 
state orbits coalesce at x=0, p — ► ±oo, making this system highly unsuitable for 
localization onto classical trajectories. We wish to examine whether another form of 
equivariant localization, such as ( [2.23j ), may apply to systems of this type. 

In considering bound states it will again be useful to change variables (of the time 
independent modes only!). We define hyperbolic coordinates 

p = | r | sinh a , 

2e 2 (4.2) 



x 



r|r|cosh 2 a 



so that H = —\r 2 . The original phase space maps into — oo < r, a < oo. The 
hamiltonian vector has only one nonvanishing component: 

X a = X = -^r 3 cosh 3 a . (4.3) 

The Lie derivative constraint Cngab = has exactly the same form as (3.4), with 
(r, 9) replaced by (r, a). Thus the general solution for has precisely the form 
(3.5). However in the present case, since a is a hyperbolic coordinate, we do not 
encounter a single- valuedness problem in defining g ra . Explicitly: 



12e 2 /(r) 



sinh a , , , . 
+ 2 tan (sinh a) 



2cosh 2 a 2 



+ AM , (4.4) 

X 
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which is perfectly well-defined. 
A trivial calculation gives 

f'(r) 

n ar = . (4.5) 

X 

The curvature can then be computed rather easily by exploiting the relation 

£2 ar . r 2_R Q , rQr ^ . (4.6) 



One obtains 



Rarar ~ ^ (f ^ f /A • (4.7) 



With these results, one can easily show that the localization formula ( |2.23|) can be 
reduced to the following simple expression: 

00 . r°° , l d ( c ii r 2 T 



— oo J — oo 



Z(T)= dr da--[——e^ rI , (4. 

Xdr V sin cl 



where the function c(r) is defined by 



cW = ■ (4 - 9) 



In fact ( J4.9| ) is a general result; it follows from ( |2.23| ) for any static potential which 
is unbounded below. 

As for the harmonic oscillator, we have obtained a remarkably simple result, which 
is unfortunately ambiguous due to its dependence on the (as yet) arbitrary function 
c(r). Actually the ambiguity appears to be much worse, since the harmonic oscillator 
result only depended on c(0). We can improve the situation by performing some addi- 
tional manipulations on ( f4.9|) . We perform the a integration, and expand sin cT as a 
power series. Then, assuming that we can interchange integrations and summations, 
we have a simple expression for the integrated partition function: 

POD °° POO 1 1 

/ dTZ{T) = -12^ Y, dr-- (4.10) 

Jo n=x J -°° r (2rc-i)-yy 

We then make the mild assumption that the only singularities of the integrand in 
( |4.10| ), assumed to be an analytic function of complex r, are poles on the real axis. 



We can then convert ( f4.10| ) into a contour integral, adding the prescription that the 



contour includes all poles except the one at r = (if we include all the poles, the 
integral simply vanishes). However, since the contour can be closed above or below, 
this is the same as including only the pole at r = 0. We therefore obtain 

r°° 00 d 3 ( 1 \ 

/ dT Z(T) = 4vrV* V lim f- r - r ^ . (4.11) 
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Thus we see that the result depends only on the behavior of c(r) as r goes to zero. 
To resolve this remaining ambiguity, we need an additional input. To see what this is, 
let us consider the main difference between the WKBJ localization of path integrals 
and the localization procedure which led to ( |2.23| ) and ( [4.91 ). In the former case 
the original path integral is reduced to a sum of small fluctuation integrals around 
the classical trajectories. In the latter case, the original path integral is reduced 
to a sum of small fluctuation integrals around individual points in the phase space 
parametrized by the time independent modes (x,p). From this fact we conclude that 
the contribution to (|4.9| ) from the region x — > oo, p ~ 0, should be approximately 
equal to the corresponding contribution in (3.39) for the free particle, provided that 
T is not near oo. 

We have applied this matching condition by computing the general metric for the 
1-d hydrogen atom in (x,p) coordinates. We then matched to the free particle metric 
( |3.14|) to leading order in large x, and to leading and next-to-leading order in small 
p 2 x. This is sufficient to give the leading order behavior of c(r) for small r: 

3 

c(r) = ^ (l + cir + c 2 r 2 + . . .) . (4.12) 



Actually, the evaluation of (|4.11|) requires knowledge of the corrections C\ and C2- 
While these can be computed in principle via perturbation theory, we have not done 
so. Regardless, the leading order behavior is enough to see that ( |4.11|) reproduces the 
correct form for the bound state spectrum of the 1-d hydrogen atom: 

e 4 

E n oc — . (4.13) 

n A 

We consider this a nontrivial success of equivariant localization techniques applied 
to physically relevant path integrals. While we are still operating at a rather primitive 
level compared to standard WKBJ, we believe that the potential for obtaining new 
results via these more general methods of equivariant localization is much greater. 
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